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Abstract 

We establish the existence, uniqueness, and stability of the stationary solution of the one-dimensional viscous 
Burgers equation with the Dirichlet boundary conditions on a finite interval. We obtain explicit formulas for 
solutions and analytically determine the Lyapunov exponents characterizing the asymptotic behavior of arbitrary 
solutions approaching the stationary one. 
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Introduction 

Burgers’ equation has the same nonlinearity form as the Navier-Stokes equations [1], It is often used as a model 
equation in studying computational methods for solving partial differential equations (PDEs) [2]. In this paper we 
establish the existence, uniqueness, and stability of the stationary solution of the one-dimensional viscous Burgers 
equation (1) on a finite interval with the Dirichlet boundary conditions (4). We use the Cole-Hopf transformation 
to give the result for any combination of A and B in the boundary conditions (4). Using a different method 
(linearization) H.-O. Kreiss and G. Rreiss (1985) gave a similar result for a subset of cases: A > |Z?|, B < 0 < A, as 
well as for Burgers’ equation with forcing [8]. We obtain explicit formulas for solutions and analytically determine 
the Lyapunov exponents characterizing the asymptotic behavior of arbitrary solutions approaching the stationary 
solution with the same boundary conditions (4). 


1 Explicit formulas for stationary solutions 

The viscous Burgers equation is the nonlinear partial differential equation 

lit + UUx = VU XX 

with v = const > 0. If we set u, to zero, for the stationary solution u = u s (x) we obtain 

uu x = vu xx . 

We note that uu x = -(u 2 )' x , therefore (2) gives 

2 vu x = ic + Cq. 


( 1 ) 

( 2 ) 

( 3 ) 


First, assume that u x 4 0 and Cq is negative, Co = - a 2 < 0 (i. e. 2vu x < u 2 ). We have dx = 2 vdu/(ir - a 2 ). 


ax , / 

— = lnlCi 

a — u 
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a + u( 0) 


a - u(0) 
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If, in addition, \u\ < a (i. e. u x < 0), then 


Ci - exp(ax/v) 
Ci + exp (ax/v) 


-2vko tanh(£o(-Y - xo)), 


where ko = 


2v’ 


XQ 


— artanh 

ko 


u( 0) 
2vk 0 ' 


while if \u\ > a (i. e. 0 < 2 vu x < ir ), then 


Ci + exp(ax/v) 
Ci - exp (ax/v) 


-2vko coth(A:o(x - xo)), 


where ko - —, xo 
2v 


— arcoth 

ko 


u( 0 ) 

2vk Q ' 


Now assume that Co is positive, Co = a 1 > 0 (i. e. 2 vu x > u 1 ). Then dx = 2vdu/(ir + a 2 ), 

cix u u( 0) , /ax \ 

— = arctan —h Ci, where Ci = — arctan-; hence u - a tan-Ci ; 

2 v a a \2v ) 

or, equivalently, 

a 1 m(0) 

u = -2vko cot(k()(x - xo)), where ko = —, xo = — arccot- 

2v ko 2vko 

Finally, if Co = 0, then u = -2 v/(x - xo); if u x = 0, then u — const and ir = |Co| = const. For convenience, all 
explicit formulas for stationary solutions are listed together in Table 1 (left column). 


Table 1. Stationary solutions u s of Burgers equation and the corresponding solutions tp s of the heat equation (6). 
H = 2v(B - A) - IAB\ 2vko = |Cq| 1/ ^ 2 , where Co is the constant in (3); u s (x ) = -2v(ln |i/> s (x, t)\)' x . 


Solution u s (x) of (1) 

Conditions on a, u x 

Conditions on A , B , H 

Solution i p s (x, t ) of (6) 

(a) -2vko cot(fco(x - xo)) 

0 < u 2 < 2vu x 

A < B, H > 0 

C sin(ko(x - xo)) exp (-vk^t) 

2vko tan(ko(x - xj)) 

(the same conditions and same solution as above) 

C cos(ko(x - xj)) exp (-vk^t) 

(b) -2v/(x-x 0 ) 

0 < u 2 = 2vu x 

A < B, H - 0 

C(x - x 0 ) 

(c) -2vko coth(ko(x - Xo)) 

0 < 2 vu x < u 2 

A<B, H < 0 

C sinh(ko(x - Xo)) exp (vk^t) 

(d) ±2vko = const 

u x = 0 

A = B 

C exp (vkgf + kox) 

(e) -2vko tanh(ko(x - xq)) 

U; < 0 

A> B 

C cosh(ko(x - xq)) exp(vkgf) 


We will now consider Burgers equation (1) with the Dirichlet boundary conditions on the interval x e [0, /]: 

u( 0, t) = A, u(l, t ) = B, (4) 

where A and B are constants. Let us find out which explicit formulas (Table 1) can represent the stationary 
solution li s of equation (1) with boundary conditions (4). Here we are concerned exclusively with solutions that 
are continuous, bounded, and sufficiently smooth everywhere on the interval x e [0, /]. 

Clearly, when A > B, the stationary solution u s can only have form (e) -2vko tanh ko(x - xq) which is the only 
decreasing function in the left column of Table 1. When A - B, the stationary solution u s can only have form (d); 
all other explicit formulas for u s defined on [0, /] are either strictly decreasing or strictly increasing functions of x. 

To examine the stationary solution u s (x) for the trickiest case, A < B, we introduce the quantity 

H = 2v(B - A) - 1AB. 

An elementary calculation shows that u s has form (b) if and only if H = 0, A < B. It remains to analyze the 
situations that yield solutions (a) and (c). We note that, at any given point (x, m(x)), any graph u s of form (a) is 
steeper than (b), while any graph of form (c) is less steep than (b). Indeed, for any stationary solution u s we have 
a constant value of Co = 2 vu x - m 2 ; solutions (a) are obtained from (3) when 2 vu x > u 2 (steeper graphs, Co > 0, 
H > 0), while solutions (c) are obtained from (3) when 2 vu x < ir (less steep graphs, Co < 0, H < 0). Thus when 
A < B and H > 0, we can only have u s given by formula (a); when A < B and H < 0 we can only have it s given 
by formula (c). 

Note also that we have not yet proved that a stationary solution satisfying boundary conditions (4) exists for an 
arbitrary combination of A and B. (We will prove this in Section 3.) Still, in the simple cases (b) A < B, H = 0 
and (d) A = B, it is already obvious that such stationary solutions do exist. 


2 






2 The Cole-Hopf transformation 


Burgers equation (1) is a rare example of a nonlinear PDE that can be linearized using a simple transformation. 
Specifically, if in equation (1) we substitute 


u(x, t ) = -2y(ln \tp(x, t)\)' x (5) 

then for the unknown function <p(x, t) we obtain the heat equation 

<Pt = vip xx . ( 6 ) 

The substitution (5) is known as the Cole-Hopf transformation [1, 2, 5, 6]. Let us discuss some interesting proper¬ 
ties of this transformation. 

Firstly, transformation (5) can produce the same solution u(x, t) of (1) from many different solutions tp(x, t ) of (6); 
these ip(x, t) may differ from each other by an arbitrary nonzero mutiplier C. Indeed, fin \ip\)' x = fin \Cip\)' x for any 
constant C + 0. 

Secondly, zero values of tp(x, t ) are mapped by (5) into discontinuities of u(x, t ). Therefore, to get a continuous 
nix, t), it is not enough to start from a continuous solution ip(x, t) of (6). We, moreover, need to restrict ourselves 
to those solutions cp(x, t ) that are nonzero everywhere on [0, /] for all t > 0. 

Further, stationary solutions u s (x ) of f 1) correspond to solutions ip s (x, t ) of (6) that may or may not be stationary. 
Explicit formulas for those < p s (x, t) that yield stationary solutions u s (x) are listed in the right column of Table 1. 
Interestingly, among these (p s (x, t) we find “non-physical” solutions of the heat equation that grow infinitely large 
when t — > oo. 


3 Existence and uniqueness of the stationary solution 

Using the Cole-Hopf transformation (5), we will now establish the existence and uniqueness of the stationary 
solution of (1), (4) for any A and B. Note that (5) transforms the problem (1), (4) into the following problem for 
heat equation (6) with the Robin boundary conditions: 


<ft = vip xx 


<Px( o, t) + ^-ipi 0, t) = 0, (fiJJ, t ) + t ) = 0. (7) 

2y 2y 

Denote by ip s the solution of (6) that under transformation (5) yields the stationary solution u s of (1). Our ip s 
must have the form cp s (x, t) = X(x) ■ T(t). (This can be checked directly by substituting ip s into (5), or simply by 
inspection of the right column in Table 1.) Here X(x) is a function of the x coordinate only, and T(t) is a function 
of time t only. Substituting this tp s into the heat equation(6) and dividing through by vTX, we get 


r _ x" 
vr ~ ~x 


= -A. 


(One ratio is a function of t only, while the other ratio is a function of x only. In order for these two ratios to be 
equal, they both must be equal to a constant which we denote -A.) 

For the function X(x), problem (6), (7) translates into an eigenvalue problem (a Sturm-Liouville problem) with 
Robin boundary conditions: 

-X"(x) = AX(x) (8) 

X'(0) + A X (0) = 0, X'(l) + ^X(l) = 0; (9) 


and for the function T ( t ) we readily obtain 


T(t) — C exp(-y/lf). 


( 10 ) 
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For it s to be continuous, p s must be nonzero everywhere on the interval [0, /]. So the question now is: how many 
eigenfunctions of (8), (9) are nonzero everywhere on [0, /]? The answer is well known: for any A and B, there is 
one and only one such eigenfunction. This follows from the familiar fact that, for any A and B in problem (8), (9), 
all eigenvalues A, (Ag < Ai < ...) have multiplicity 1, and the respective eigenfunction X,(x) has exactly i zeros 
inside the interval (0, /); see [3, pp. 14-18]. Thus, in problem (8), (9) we are interested in the eigenfunction Xg(x) 
that has no zeros for x e [0, /] and corresponds to the least eigenvalue do. For p s we find, up to a nonzero multiplier 


C, 


i p s (x, t ) = CXq(x) ■ exp(-vdo0 (^ S has no zeros for x e [0, /]). 


Therefore, for any A and B, there exists a unique stationary solution u s (x) of Burgers equation (1) with boundary 
conditions (4): 


u s (x) = -2v(ln \p s {x, t)\)' x = -2v(ln|X 0 (x)|)' t . 


4 Stability and Lyapunov exponents 


Now let us study the evolution of the absolute value \u - u s | for an arbitrary non-stationary solution 

u(x,t) = -2v(ln \p{x, 01)» = -2v ^| , 

(p(x, t) 

where both u(x, t) and u s (x) satisfy the Burgers equation (1) with boundary conditions (4), and < p(x, t) is a suitable 
positive solution of (6). It is known that the solution u(x, t ) exists for “reasonable” combinations of the boundary 
conditions (4) and initial condition u(x, 0) [9]. We say that u s is stable if \u - u s | —> 0 as 1 —> oo, for an arbitrary u 
obeying (1), (4). We have 


|u - ir | = 2v 


= 2v 


= 2v 


fx_ _ f* 

ps p 

<fi(tPx - Vx) + <Px(<p - W S ) 


W S x -‘fS'Px 


ipS (p 


P)S ip 


= 2v 


<PxV - Wx 


ips ip 


Here we have introduced the notation ip — ip - tp s . Taking into account that u = -2vp x /tp, for all x e [0, /] and all 
t > 0 we obtain the estimate 



+ 2v 

$x 

ipS 


ipS 


< max |w(x, 0)1- 


+ 2v 

$x 

VG[0,/] 

(pr 


ipS 


( 11 ) 


In inequality (11) we have used the maximum principle for Burgers equation: the solution u(x, t) attains its maxi¬ 
mum either in the initial value u(x, 0) or at the boundary of the interval [0, /]. (A discussion of maximum principles 
for PDEs can be found in [4, 7, 9]. The proof of the maximum principle for Burgers equation is similar to that for 
linear parabolic PDEs.) 

Expand p(x, t) in a series over the system of eigenfunctions Xj(x) of (8), (9): 


OO oo 

tp(x, t) = Yj a iXi(x)Tj(t) = ^ Pi(x, t ), Tj(t) = exp(-V/l,f), T,{ 0) = 1. (12) 

i =0 i—0 


In this series, the term po(x, t) = agXg(x)Tg(f) is the same as p s (Table 1) up to a constant nonzero multiplier. Let 
us choose C in the expression of p s (Table 1) so that pg - p s . If we now compute the difference p - p s , the term 
po(x, t ) will cancel out, and we get 

OO 

<t> = <P - <p s = 2 <Pi(x, t). (13) 

i= 1 

Since T,(t) = exp(-V/l,f). we see that p\(x,t) becomes the largest term in (13) when t —> oo (assuming a i 4= 0 in 
(12)). We then have 

max x expI-vdoO, max \ip\ x exp(-vdif), max \p x \ x exp(-vdif) as f —» oo, 

-v e [0,/] x e [0,/] xe[0,/] 
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so the estimate (11) results in 


max | u - u s \ x exp(-v(/li - Ao)t) as t —> oo. (14) 

x e [0,/] 

This paves the way to proving the stability of the stationary solution u s . Indeed, the difference | u - u s \ is an 
exponentially vanishing quantity when t —* oo. Nevertheless, the convergence of \u - u s \ to zero might turn out to 
be very slow; this is the case when the least two eigenvalues To and A\ in problem (8), (9) differ only slightly. 

We got the estimate (14) under the assumption that a i 4 0 in (12), that is, in the series expansion of p over the 
system of eigenfunctions Xfx) there is a nonzero term ip\ containing the eigenfunction X\ (x). However, if it so 
happens that one or more initial terms in (13) are zero, then the series (13) for ip = <p> - <p s will start at some ip„ 
(n > 1). In the general case, therefore, instead of (14) we would have 

max | u - u s | - exp(-v(/l„ - Ao)t) as t —> oo, (15) 

x e [0,/] 


where n is the number of the first nonzero term in the series expansion of ip — ip - ip s (13). We have thus proved 
that the stationary solution u s is stable: \u - u s \ —> 0 as t —> oo. 

Note that the functions ip, (i = 1,2,...) in (12) have the same explicit formulas as ip s (Table 1), except that each <pt 
contains its own values in place of ko and xq; let us denote these new constant values by ki and Xj, respectively. 

All constants ki and x, can be found if we substitute the general solutions of (8) (trigonometric, exponential or 
hyperbolic functions) for the eigenfunctions Xfx) ( i = 0,1,2,...) in the boundary conditions (9). In most cases 
(i. e., cases (a), (c), (e) in Table 1), this substitution yields the following transcendent equations for£; = kf: 


■- j + qfi for pi of forn 

hi 

-- j - qfi for pt of forn 
hi 

where = kf > 0, p 


— kf > 0, and 

(16) 

, A, = -kf < 0, 

(17) 


1AB 


2v 


q ■ 


2 v(B-A) 1 l(B-A) 


The transcendent equation (16), with cot(-',, may correspond to any i, whereas equation (17), with coth (■), may 
correspond only to i = 0,1 (the least two eigenvalues To, Ti) because hyperbolic functions cannot have more than 
one zero value on the interval [0, /]. 


When A = B in (4) and (9), we have an exceptional case: all kj and T, can be found in a closed form. Here the 
interval [0, /] contains a whole number of semiperiods of the eigenfunction Xfx) = sin (kfx - x,)), i = 1,2,..., 
which readily yields 


kj = 


ni 

7 


a = i,2,...), 


while 



Xo(x) = C exp(+fco.r); see Table 1 (d). 


Therefore, if A = B, we find 


A 


n 



(n > 1 ), 



and A n - To 



cf. (14), (15). 


Now we will reuse the customary definition of Lyapunov exponents in the context of problem (1), (4) for Burgers 
equation. Let u(x, t) be a solution of (1),(4). The Lyapunov exponent pi of this solution is defined as 


p = lim sup 

t —>oo 



(18) 


This definition, in general, depends on our choice of the norm || • ||. If u(x, t ) behaves so that \\u - u s || x exptcif) as 
t —> oo, then it is easy to see that 6 is the Lyapunov exponent of this u(x, t). 

Let us use the norm defined as the maximum absolute value: 


||w(x)|| = max |w(x)|. 

X £ [0,/] 
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Then estimates (14), (15) allow us to determine all Lyapunov exponents for any u(x, t) satisfying (1), (4): 

Hi = -v(Aj - To), i= 1,2,..., (19) 

where, as before. A, are eigenvalues of (8), (9). Solutions u(x, 1) corresponding to the Lyapunov exponents p, can 
be written simply as 

Uj(x, t ) = -2v(ln | tp s (x, t ) + tpAx, t)\)' x , i = 1,2,..., 
where ipfx, t) is the respective term of (12). For example, when u s has the form (a) in Table 1, we have 

ip s (x, t) — C sin(ko(x - xo)) exp(-v^/) (<p s has no zeros for x e [0, /]), 

ipi(x, t) - ai sin(k,(x - x,)) exp (-vkft) (ip, has i zeros for x e [0, /]), 

and we can write a solution ufx, t) corresponding to the Lyapunov exponent /./, as follows: 


Ui(x,t ) = —2v 


Cko cos(£o(x - xo)) + ajkj cos(k,(x - x,)) • exp(-v(kr - kfyt) 
C sin(^o(x - xo)) + a,- sin(fe,(x - x,)) • exp(-v(^ - k^)t) 


( 20 ) 


Because each individual term in series (12) satisfies the Robin boundary conditions (7), each function m,(x, t) 
defined as above must satisfy the Dirichlet boundary conditions (4). 

We have thus determined the Lyapunov exponents in the nonlinear problem (1), (4) for Burgers equation: we have 
found that formula (19) relates the Lyapunov exponents /.(, to the eigenvalues A, of the linear problem (8), (9). 
All Lyapunov exponents p, are negative; there are countably many of them; we can write explicit formulas for 
the corresponding solutions u,(x,t) of Burgers equation (1). This is an interesting example of a situation where 
one can analytically determine the Lyapunov exponents for solutions of a nonlinear PDE with Dirichlet boundary 
conditions. 
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